Coherence of nitrogen-vacancy electronic spin ensembles in diamond 
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We present an experimental and theoretical study of electronic spin decoherence in ensembles 
of nitrogen-vacancy (NV) color centers in bulk high-purity diamond at room temperature. Under 
appropriate conditions, we find ensemble NV spin coherence times (T2) comparable to that of single 
NVs, with T2 > 600 /is for a sample with natural abundance of ^^C and paramagnetic impurity 
density 10^^ cm~^. We also observe a sharp decrease of the coherence time with misalignment of 
the static magnetic field relative to the NV electronic spin axis, consistent with theoretical modeling 
of NV coupling to a ^^C nuclear spin bath. The long coherence times and increased signal-to-noise 
provided by room-temperature NV ensembles will aid many applications of NV centers in precision 
magnetometry and quantum information. 



The Nitrogen- Vacancy (NV) color center in diamond 
is a promising solid state platform for studying the quan- 
tum dynamics of spin systems. Experiments demonstrat- 
ing spin-state control and manipulation using single NV 
centers and their nuclear environmenlP'^ have inspired 
recent work towards quantum information applications 
of NV diamond, as well as applications to precision mag- 
netic field sensing,^ ^ including the life sciences.^ Com- 
mon to all proposals is a desire for long NV electronic 
spin coherence times, ideally at room temperature, and 
high signal-to-noise of the NV spin-state-dependent op- 
tical fluorescence. 

Electronic spin decoherence of the NV center (S* = 1) 
is governed by interactions with the surrounding bath of 
nuclear and paramagnetic spins (Fig. |l(a)] ). In high pu- 
rity diamond, decoherence is dominated by ^^C nuclear 
spins (/ = 1/2),^^ which are dispersed through the crys- 
tal with a natural abundance of 1.1%. Electronic spin 
coherence times (T2) longer than 600 /is at room temper- 
ature have been observed for individual NV centers in 
this type of samplc^^^ T2 can be increased by isotopi- 
cally enriching the diamond^ for example, T2 > 1.8 ms 
at room temperature has been observed in isotopically 
enriched ultra-pure diamond with 0.3% ^^C.^ 

In addition to long electronic spin T2, enhancing the 
signal-to- noise of NV spin-state-dependent fluorescence is 
of significant benefit to precision magnetometry. To th at 
end, ensembles of NV centers have been proposed,^^^ for 
which signal-to-noise ideally increases as the square root 
of the number of NVs in the ensemble. The use of large 
spin ensembles with good coherence times is also impor- 
tant for collective quantum memories .'^^'^I^ Increasing the 
density of NV centers, however, is accompanied by an 
increased density of residual nitrogen paramagnetic im- 
purities, which can become the dominant source of NV 
spin decoherence, reducing T2 below the limit set by ^^C 
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FIG. 1: 



(a)| NV electronic spin axis is defined by nitrogen and 

NV 



vacancy sites, in one of four crystallographic directions, 
orientation subsets in an ensemble can be spectrally selected 
by applying a static magnetic field, Bq. Also shown are ^^C 
nuclear spins and othe r pa ramagnetic impurities such as sub- 
stitutional nitrogen. (b) NV center electronic energy level 



structure. Spin polarization and readout is performed by op- 
tical excitation and fluorescence detection. Ground state spin 
manipulation is achieved by resonant microwave excitation. 
The ground sta te t riplet has a zero magnetic field splitting 
D :^ 2.87 GHz. (c) Hahn-echo experimental sequence used in 



present measurements. 



nuclear spins alone Furthermore, the large variation 
of coherence times for individual NV electronic spins, due 
to the random distribution of ^^C near each NV center, 
can also affect ensemble coherence. Therefore, a detailed 
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understanding of decoherence mechanisms affecting NV 
centers - both single NVs and ensembles - accompanied 
by the development of control techniques to mitigate de- 
coherence is of great importance for both precision mea- 
surement and quantum information applications. 

Only limited studies of electronic spin decoherence in 
ensembles of NV centers have been previously reported. 
At room temperature, ensemble T2 times of 58 jas were 
measured in CVD diamondpl whereas at low temper- 
ature (< 2K) and high magnetic field (> 8T), NV en- 
semble T2 of 250 /is was measured in high-temperature 
high-pressure type lb diamond, decreasing to T2 < 10/is 
for temperatures > 20 K.^ 

In the present Rapid communication, we report an ex- 
perimental and theoretical study of the coherence proper- 
ties of NV electronic spin ensembles in room temperature 
diamond samples of different paramagnetic nitrogen (and 
consequently NV) concentrations. For a lower nitrogen 
density sample (~ lO^^cm"^) with natural abundance of 
-^^C, we find NV ensemble T2 in excess of 600/is, compara- 
ble to the best results for single NV center measurements 
in natural isotopic abundance diamond and an order of 
magnitude greater than previous room-temperature en- 
semble measurements. For a higher nitrogen density sam- 
ple 5 X 10^^ cm~^) with natural abundance of ^^C, we 
find an NV ensemble T2 ~ 300 jas. Furthermore, for both 
samples we find a sharp decrease in the NV ensemble T2 
with misalignment of the static magnetic field relative to 
the NV electronic spin axis being studied, which is con- 
sistent with our theoretical modeling of an ensemble of 
NV electronic spins interacting with a ^^C spin bath.l^ 

For high-purity diamond, decoherence of single NV 
electronic spins is dominated by hyperfine interactions 
with nearby ^^C nuclear spins. The contact term of this 
interaction decreases exponentially with separation; af- 
ter a few lattice sites it is not larger than a few MHz.^^ 
Meanwhile, the dipolar part of this interaction decreases 
as and is responsible for the collapses and revivals 
observed in Hahn-echo measurements of single NV cen- 
ters .^^ When the externally applied static magnetic field is 
aligned with the NV axis, the dipolar field contributions 
of all ^^C nuclei cancel after each 27r Larmor precession 
period of the ^^C nuclear spins, as can be measured by the 
Hahn-echo sequence. Over longer timescales (> 600 /is 
for natural abundance ^^C), weak dipole-dipole interac- 
tions between ^^C nuclei and between ^^C and paramag- 
netic impurities induce fluctuations ("flip-flops") in the 
NV-^^C hyperflne interaction, leading to NV spin de- 
coherence. The predicted form of this decoherence for 
the Hahn-echo signal of an individual NV is exp [— /3r^], 
where n is between 3 and 4 and the constant /3 depends 
on details of the relative location of nearby ^^C nuclei.^ 
However, when the magnetic fleld makes a small angle 
with the NV axis, the Larmor precession frequency Qi of 
individual ^^C nuclei is modifled by the hyperflne interac- 
tion with the NV center to become dependent on and 
the relative NV-^^C position r^: Qi = -\- 5Q{0^ri). 
As a consequence, ^^C nuclei do not all precess with the 



same frequency and thus their modulations of the NV 
spin coherence do not rephase at the same time in the 
Hahn-echo sequence, inducing NV spin decoherence of 
the form exp [—a (0) 6r^] .^^ Here a (6) describes the mis- 
alignment angle dependence of the imperfect Hahn-echo 
rephasing, and b is proportional to the square of the hy- 
perflne interaction between the NV electronic spin and 
the nearest ^^C nuclear spin. For small 6, a{0) 0'^^ 

Note that the physics of this angle-dependent NV spin 
decoherence is fundamentally different from previously 
observed effects for do nor electron spins in semiconduc- 
tors such as SiiP.^J^For Si:P and similar 5 = 1/2 sys- 
tems, the electronic spin quantization axis points along 
the external magnetic fleld; whereas for NV centers, 
the anisotropy of the electronic distribution causes spin 
quantization along the NV axis for small magnetic flelds 
(< 500G). In Si:P and similar systems, electron spin 
decoherence is dominated by fluctuations of a dipolar- 
coupled nuclear spin bath, and is maximized when the 
magnetic fleld is aligned along the [111] axis, due to the 
enhanced flip-flop rate of the nuclear spins arising from 
the angular dependence of their dipolar coupling. As 
noted above, NV spin decoherence is minimized when 
the magnetic fleld is aligned with the NV axis. 

In diamond samples where the concentration of NV 
centers is low enough that they do not interact signifl- 
cantly with each other and with nitrogen paramagnetic 
impurities, the ensemble Hahn-echo signal can be con- 
sidered as the average of many independent signals from 
individual NV centers. In this case, the ensemble Hahn- 
echo signal envelope, E{t)^ can be signiflcantly different 
from measurements of single NV centers, because an 
ensemble contains a broad distribution of spin coherence 
lifetimes due to variations in the location of ^^C nuclei 
proximal to individual NVs in the ensemble: 



E(t) = J dbexp {-a {0) br^) f 



(6). 



(1) 



Here, / (b) is the probability distribution for the mag- 
nitude of the NV-proximal ^^C hyperflne-interaction- 
squared in an ensemble. Depending on details of this 
distribution, the ensemble Hahn-echo signal envelope can 
have either a sing le exponential or Gaussian decay at 
large times r P^ ^ ^ ^ 

In our experiments we used a custom-built, wide-fleld- 
of-view fluorescence microscope to measure the coher- 
ence properties of large ensembles of NV centers. The 
electronic spin state of negatively charged NV centers 
can be polarized, manipulated and rea d-out using optical 
and microwave excitation (see Figure fT^bJ). Spin-state- 



dependent radiative-relaxation enables NV polarization 
and readout: NVs initially in = optically cycle be- 
tween the ground and an electronically excited "bright" 
state, while NVs initially in \ms\ = 1 have a signiflcant 
branching to a long lived metastable "dark" state that 
non-radiatively relaxes to = 0. Optical excitation is 
provided by a switched 532 nm laser source, focussed onto 
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FIG. 2: Example Hah n-echo measurements of ensemble NV electronic spin decoherence and depend ence on static magnetic 



field orientation: 
density ^ 1 x 10^ 



sample A with NV density ~ 



X 10^"^ cm""^ and N density ^ 10^^ cm""^; and |(b)| sample B with NV 
and N density ^ 5 x 10^^ cm""^. Clearly seen in this data is a sharp decrease of the ensemble NV 
coherence lifetime T2 as the Bq field is misaligned from the NV electronic spin axis, and as the paramagnetic impurity level 
increases. For each data plot, the left vertical axis indicates the measured probability for NVs in the ensemble to be in the 
nis = state at the end of a Hahn-echo experiment, as a function of the echo time r; the black line is a fit of exp [— (t/t2)^] 
to the overall Hahn-echo signal envelope, where the exponent n is a fit parameter and the long term baseline is taken to be 
0.5 (the mixed spin state); and the errors in 0, T2 and n are given by one-standard-deviation confidence intervals of the fit. 
The red line is a fit to the observed echo-signal modulations, which are collapses and revivals of NV spin coherence induced by 
precession of the natural abundance ^^C nuclear spin bath in applied, hyperfine and dipolar magnetic fields. 



the sample with a 20x objective. Fluorescence emanating 
from the sample is collected back through the same ob- 
jective, then filtered and imaged onto a CCD array. The 
use of a CCD array allows the sample response to be spa- 
tially resolved with resolution < 1 /im and field-of-view 
> 100 /im. The ensemble of NVs is coherently manipu- 
lated on a ground-state electronic spin transition with a 
resonant microwave magnetic field generated from a loop 
antenna placed close to the diamond sample, which pro- 
duces a homogenous microwave Bi field over the region 
of interest. The degeneracy of the = to = ±1 
ground-state spin transitions is lifted by a uniform static 
^0 magnetic field, applied using a permanent magnet 
trimmed by a 3-axis set of Helmholtz coils. Precise align- 
ment of ^0 allows a subset of NVs in the ensemble, cor- 
responding to one of four crystallographic orientations, 
to be spectrally distinguished by their spin transitions. 
The orientation and magnitude Bq of the applied mag- 
netic field relative to the axis of one class of NV centers 
is determined from the Zeeman shifts of the NV ground- 
state electronic spin resonances (ESR) together with the 
measured Hahn-echo revival frequency, Bq = ^jl/Ji^c^ 

We employed two CVD diamond samples in our mea- 
surements of NV ensembles. Each sample consists of a 



pure diamond substrate (3 x 3 x 0.2 mm) with < 1 ppm 
of ambient nitrogen impurity and < 10^^ cm~^ NV den- 
sity, which was overgrown with a high NV content layer 
~ 2 /im thick. The density of NV centers in the thin layer 
region of sample A and B is approximately 3 x 10^^ cm~^ 
and 1 X 10^^ cm~^ respectively, as determined by fluores- 
cence measurements normalized to single NV measure- 
ments made on each sample using a confocal microscope 
optimized for NV measurements, similar to that used in 
Ref. [71 The density of nitrogen impurities in the samples 
is approximately 50 times larger, as determined by the 
CVD growth procedure. A Hahn-echo pulse sequence (il- 
lustrated in Fig. |l(c)[ ) was used to measure the electronic 
spin coherence time of the NV ensemble in each of the two 
samples. NV centers were first prepared in the = 
ground state by applying a 532 nm laser excitation pulse. 
The Hahn-echo sequence {71/2— r — 7r — r — 71/2) was then 
repeated for a range of free spin evolution periods r (the 
final 71/2 pulse in this sequence projects a coherent super- 
position of NV spin states into a population difference, 
which is detected via spin-state-dependent fluorescence) . 
The data presented here is the integrated response of re- 
gions with area ~ 30/im^ and 2/im depth in each sample, 
averaged over ^ 10^ identical experiments. Thus, NV en- 
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sembles consisted of ~ 10^ NV spins of a single crystallo- 
graphic class in the detection region. Over macroscopic 
distances (~ 1 mm), each diamond sample was found to 
have uniform spin coherence properties. 
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FIG. 3: Summary of Hahn-echo measurements as a function 
of misalignment angle bet ween the static magnetic field, Bp, 
and the NV electronic axis: (a) coherence lifetime T2 and (b) 
decay exponent n, determined from fits of the signal enve- 
lope to exp [— (r/T2)^]- Sample A (lower NV and nitrogen 
concentration) was measured for ranging from 90 to 180 
degrees; sample B (higher NV and nitrogen concentration) 
was measured for from to 180 degrees. Solid curve in (a) 
is the prediction of the NV ensemble decoherence model, as 
described in the main text. 



As shown in Figs. |2] and [3j we measured ensemble NV 
coherence times comparable to that for single NVs in each 
sample: T2 > 600 jas for sample A and T2 ~ 300 /is for 
sample B, with the static magnetic field Bq well aligned 
with the NV electronic spin axis. Our measurements also 
confirm the predicted^ ^ sharp decrease of NV T2 with 
small misalignment of Bq from the NV axis. Example 
Hahn-echo measurements for three representative mis- 
alignment angles are shown in Fig. |2] for samples A 
and B; whereas Fig. [3] summarizes the results of a larger 
set of systematic measurements of ensemble NV T2 as a 



function of ^. In all these measurements, the Hahn-echo 
signal envelope was fit to ex p [— (t /t^)^] to determine 
the coherence lifetime T2 (Fig. |3(a)| ) and decay exponent 
n (Fig. |3(b)[ ). Good agreement is found between our 
experimentally-determined values for T2 as a function of 
and theoretical predictions (see solid curve in Fig. |3(a)[ ), 
which include the decoherence contributions from proxi- 
mal ^^C nuclear spins via the form exp [—a (0) 6r^] IHEni 
Decoherence due to ^^C "flip-flops" was included phe- 
nomenologically by multiplying the simulated NV echo 



signal by exp 



The low values for the 



(^/t2(0=o)) 

decay exponent n in ensemble experiments is explained 
by the large dispersion of the dominant NV-^^C dipole- 
dipole interaction (see Eq. [T]). The distribution of this 
interaction squared, / (6), was determined by calculating 
the effect of ^^C nuclei placed randomly in the diamond 
lattice at natural isotopic abundance. For the dominant 
NV-^^C dipole-dipole interaction, we find that this dis- 
tribution scales as / (b) ~ for ^^C nuclei within 
about 10 nm of an NV center. Thus, / (6) more closely 
resembles a Lorentzian than a Gaussian distribution, re- 
sulting in Hahn-echo signal envelopes that decay more 
like a single exponential than a Gaussian at large times 
(see Fig. [2]). Note that for sample B, the experimentally- 
determined decay exponent is n 1 even for ^ 
(Fig. 3(b) ), which could result from interactions between 
NV spins and the higher density of paramagnetic impu- 
rities (nitrogen and NVs) in this sample. 

In summary, we demonstrated experimentally that 
large ensembles of NV centers in high-purity diamond 
with natural abundance of ^^C can have electronic spin 
coherence lifetimes at room temperature that are com- 
parable to the best measured for single NV centers 
(T2 > 600 /is). We also found a sharp decrease in NV 
T2 as the applied magnetic field is misaligned from the 
NV axis, consistent with the predictions of our model 
for an ensemble of NV electronic spins {S = 1) coupled 
via position-dependent hyperfine interactions to a ^^C 
nuclear spin bath, leading to imperfect ^^C spin-echo re- 
vivals and hence NV decoherence. Our results demon- 
strate the potential of NV ensembles for applications in 
precision magnetometry in both the physical and life sci- 
ences, combining long electronic spin coherence times at 
room temperature with the enhanced signal-to-noise ra- 
tio provided by many NVs in the detection volume.^- In 
addition, the demonstrated techniques could be used to 
increase the coherence time of solid-state quantum mem- 
ories, such as an NV electronic spin ensemble coupled to 
a superconducting resonatoi^^ or flux qubit.^^ 
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